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p ,! Abstract 

D ■ 

,-^ I The renormalization factor and 0{a) improvement coefficient of four-quark operators are cal- 

culated perturbatively for the improved Wilson fermion action with clover term and the Iwasaki 
gauge action. With an application to the K — )■ vrvr decay amplitude in mind, the calculation is 

o, 

Tij- , restricted to the parity odd operator, for which the operators are multiplicatively renormalized 



o 



X 



without any mixing to operators that have different chiral structures. 



I. INTRODUCTION 

Calculation of hadron matrix elements of phenomenological interest is one of major appli- 
cations of lattice QCD. When one tries to calculate weak matrix elements of four quark oper- 
ators using the Wilson fermion action, however, one encounters an obstacle since unwanted 
mixings with operators having wrong chirality, which are prohibited in the continuum, is 
generally introduced through quantum corrections. 

This problem is absent for parity odd operators. Using the discrete symmetries of parity, 
charge conjugation and flavor exchanging transformations, it was shown [1] that the parity 
odd four quark operator has no extra mixing with wrong chirality operators even without 
chiral symmetry. This is a welcome feature for calculation of the K ^ Tin decay amplitude 
with the Wilson fermion action. 

An improvement with the clover term is indispensable for the Wilson fermion action. 
Since the renormalization group improved gauge action of Iwasaki type j^] has good scaling 
property already at lattice spacings around a~^ ~ 2 GeV, the combination of the Iwasaki 
gauge action and the improved Wilson fermion action with clover term is a plausible choice 
for numerical simulations. Unfortunately renormalization factors of the four quark opera- 
tors are not available for this combination of the actions. In this paper we calculate the 
renormalization factor of four quark operators which contribute to the K -^ titt decay to 
one-loop order in perturbation theory. 

This paper is organized as follows. In Sec. [Ill we briefly introduce the action and the 
Feynman rules relevant for the present calculation. In Sec. IIIII the AS = 1 four quark 
operators are introduced, which contribute to K ^ nn decay. The one loop contributions 
are briefly reviewed in Sec. II VI for gluon exchange diagrams and are calculated in Sec. IVJfor 
penguin diagrams. The renormalization factors are evaluated in Sec. I VI I for the MS scheme. 
Sec. IVIII is devoted for an evaluation of 0{a) effect. Our conclusion is in Sec. IVIIII 

The physical quantities are expressed in lattice units and the lattice spacing a is sup- 
pressed unless necessary. We take SU(A^) gauge group with the gauge coupling g and the 

second Casimir Cp = , while iV = 3 is specified in the numerical calculations. 

27V ' ^ 



II. ACTION AND FEYNMAN RULES 
We adopt the Iwasaki gauge action 

'S'giuon = — S Co ^ TrUpi + Ci ^ TrUrtg > (11- 1) 

I. plaquette rectangle ) 

with Ci = —0.331 |2| and the improved Wilson fermion action with the clover term 

•S'fermion = ^i){n) (^7/.^/. - 2^^ + "^oj V'(^) " CSW ^ ^ «fi'-^„0-^,.P^i.(n)^„,(II.2) 
n n fj,,u 

where a^,^ = | [7^,7i/]. We shall set r = 1 in the following. 

Weak coupling perturbation theory is developed by expanding the action in terms of 
gauge coupling. The gluon propagator of the Iwasaki action is given as an inverse of the 
action kernel 



Gi^ip) = {D-\p))^J^^, (II.3) 

D^uip) = Pt^Pu + ^{Pp^i^u - PpSp^)q^pPp, (II.5) 

p 

qpu = il-5^u){l-ciipl + pl)), (II.6) 

P, = 2sin^, (II. 7) 

where we adopted the Feynman gauge. 

The quark propagator is given by that of the ordinary Wilson fermion action 

^^^^^ - p^ + w{py — ' ^"-^^ 

Pp = sinp^, (II.9) 

VI/(p) = M + r^(l-cosp^). (11.10) 

We only need single gluon vertexes for our present calculation, one from the standard 
gauge coupling 

Vi'l{k,p) = -igT^ U^cos-{-kp+ p^,) - irsiYi-{-k^+ pM (H.ll) 
and another from the clover term 

yip^{k,p) = -gT^csw^ I Yl ^'^'^ ^^"^(^^ + ^^) ) ^°^ 2^P^' + ^/^)' (11.12) 



where k and p represent incoming momentum into the vertex as is shown in Fig. [H T^ 
(y4 = 1, . . . , A^^ — 1) is a generator of color SU(A^). 

III. FOUR-QUARK OPERATORS 

We shall evaluate the renormalization factor of the following ten operators 
Q(i) = (sd)^ (uu)^ , g(2) ^ (^ X ^)^ (^ X ^)^ ^ (III 1) 

Q^^^ = {sd)L{uu + dd + ss)^, Q^^^ = {s x d)^ {u x u + d x d + 1 x s) ^, (111.2) 

Q^^^ = {sd)^{uu + dd + ss)^, Q^'^^ = {-3 x d)^ (u x u + d x d + s x s) ^, (III.3) 

Q^'^^ = (sd)^ i uu dd ss j , Q^®^ = (s x d)^ iuxu dx d sxsj , 

(III.4) 
Q^^) = {sd)j^ f uu dd Is 1 , (5*-^°'' = (s X d)^ Iuxu d x d sxs) , 

(in.5) 

where 

(srf)^/^ = s7^(l±75)t/ (III.6) 

and X means the following contraction of the color indices 

Q(2) = (s X d)^ {U X U)^ = {^adbJL (UbUa) l ■ (ULT) 

We notice that these operators are not all independent, satisfying the relations 

QW = _Qa)+g(2)+g(3)^ (ni.8) 

Q''^ = \ (3Q(^) - Q^'^) , (III.9) 

g(io) = 1 (3Q(2) - qW) . (III.IO) 

We are interested in the parity odd operators, which contribute to the K ^ titi decay 
amplitude 

^(2n-l) _ _^(2n-l) _ ^(2n-l) ^(2n) _ _^(2n) _ ^(2n) /^ _ i 9 p:VTTT 1 1 "I 

'^VA+AV ~ ^VA ^AV ' ^VA+AV ~ ^VA ^ AV 1 V'' ~ ^1 ^I'-'lK^^^-^^J 

^{2n-l) _ ^(2n-l) _ ^(2n-l) ^(2n) _ ^(2n) _ ^(2n) / o 4\ /ttt ioN 

Qi.'r'^ = (^rf)v E <n^g)A, (ni.13) 

q=u,d,s 



Q^Iv-'^ = isd)^ J2 c.tHm)v. (111.14) 

q=u,d,s 

Q%i = (^ X d\ Yl "J'^^ (^ ^ ^)a ' (m-15) 

q=u,d,s 

qJ? = (3 X rf)^ ^ «(") (g X g)^ , (III.16) 

q=u,d,s 

where the coefficients are given for q = u,d,s 

a« = (1,0,0), (III.17) 

af ) = af ) = (1, 1, 1) , (111.18) 

4^^ = «f) = (l,-i,-0 (III.19) 

and current- current vertex means 

{-sd)y (qq)^ = (s7^rf) (gT^Ts?) • (III.20) 

There are two kinds of one loop corrections to these operators. One is given by gluon 
exchanging diagrams in Fig. [2] and the other by the penguin diagrams in Fig. |3l The 
gluon exchanging diagrams preserve the flavor structure and there occur no mixing between 
operators with different n. On the other hand the penguin diagrams mix any operator to 
the penguin operator unless it belongs to a different representation of the flavor SU{3)f. 

There are a number of parity odd operators having different chirality 

q=u,d,s q=u,d,s 

Of^i^s = (^ X d)s E "S"^ (^ X ^)p ± (^ X ^)p E "S"^ (^ X ^)s ' (111-22) 

q=u,d,s q=u,d,s 

0^;~-'^ = m^ E «S"^ (^^)t ' (III-23) 

Oj!;^ = (^ X d)^ E "f (^ ^ ^)t ' (III-24) 

q=M,d,s 

{sd)s {qq)p = i'sd) (575?) , (sd)^ {qq)f = (m^^d) {qa^^'^^q) (III.25) 

Chiral symmetry does not prohibit mixings of these operators in the Wilson fermion system. 
However, it was shown in Ref. l| that the Wilson fermion system has sufficient set of discrete 
symmetries to protect the parity odd operators Qva±av from such mixings. 



For the gluon exchanging diagrams, the operator mixing can be studied with the following 
operator having four flavors 

0S,±r'r = (^iV^2)r (^^3^4)^ ± (^^1^2)^, (V^3V^4)r , (ni.26) 

0^4'±r'r = (^1 X V'2)p (V's X ^4)^/ ± {i^i x ^^2)^., (V^3 x^i)^- (111-27) 

It was proved in Ref. [l| that parity, charge conjugation and two types of flavor exchanging 
transformations 

S' = (^1 O ^2, ^3 ^ ^4) , S" = (^1 ^ ^4, V^2 ^ ^3) (111.28) 

are sufficient to show that the mixings occur only between O^j^^^y and O^j^^^y, or between 
Oyj^_^y and Oy^_^y. Wc uotlce that the Fierz rearrangement leads to relations between 
the Fierz partners 

/-)(<^) — 9n(°)^ n(°) — on^'^'i^ riTT 9q\ 

^VA-AV ~ '^^SP-PS^ ^VA-AV ~ '^^SP-PS^ 1,111. zyj 

4p-P5 = (^1^^4)5 {^s^2)p - (V^iV^4)p (^^3^^2)5 . (111.30) 

Thus OvA-Av and Osp-ps are basically the same operator rearranged with each other, and 
if we include the Fierz partner Oy^_^y it also mixes with Osp-ps- 

For the penguin diagram we need to keep the three flavors structure. Hence flavor ex- 
change of 

S' =[d^ s.l^i) (III.31) 

is a symmetry. Together with parity and charge conjugation one can show that the mixing 
occurs only among Qva±av'^- 

We note that these discrete symmetries still allow mixings with lower dimensional parity 
odd operators 

{nid - rris) (s-f^d) , (III. 32) 

{nid - rris) d^ (s-ff^-f^d) , (III.33) 

{rrid - rris) i^F^^a^y-i^d) , (III. 34) 

(mrf - nis) ysF^yO^yd) (III.35) 

proportional to a mass difference (m^ — rris). 



IV. ONE LOOP CORRECTION FROM GLUON EXCHANGING DIAGRAMS 

We consider the following four fermi operator 

O^y = {T^'^)^,.^^, {Tx ® Ty)^^.^^, {sa,.d,,^) {q,^^q',^s) , (IV.l) 

in order to evaluate the one loop correction from the gluon exchanging diagrams given in 
Fig. [21 where T represents the color factor 

(^^'"-'l,e. = (1^1).,.. = ^^^^c, (for Q(--)), (IV.2) 

(^^"l,e. = (101).,.. = "<= (fo^ ^^"^)- (IV.3) 

and r is the gamma matrix 

(rx ® ry)„,^^, = (r^),, (r^)^, , (iv.4) 

Ty = 7^, Ta = 7m75 (IV. 5) 

where summation over /i is taken. ^ 

Since the one loop correction has already been evaluated in Ref. |3i] for various gauge 
actions, we just briefly review the result for the Iwasaki gauge action. We consider one loop 
corrections to the amputated four quark vertex 

Ik;XY = (O^XYSaJ-bpqc^q'dSJ ■ (IV.6) 

The contributions from the diagrams (a) and (a') are given by 

4:iy = Jk^ f (^^im(0, 1)Sf{1)TxSf{1)V,.{-1, 0) ® TyG,,{l), (IV.7) 

45y = Jt^'^x ® f_ ^^i.(0> l)Sp{l)TySp{l)VU-l, 0)G,M)- (IV.8) 

The contributions from the diagrams {b) and (&') are 

4?^ = Jk'^ f ^,v,,{o,i)SF{i)rx®ryS^{i)vu-i,o)G,^{i), (iv.9) 
^^^y = Jk^ f ^rx5^(/)Vi,(-/,o) ® vuo,i)SF{i)ryG,,{i). (iv.io) 

The contributions from the diagrams (c) and (c') are 

4:L = Jt^ f_ (^^im(0, 1)Sf{1)Tx ® V^i.(0, -l)Sp{-l)TyG,Al), (IV.ll) 

43y = Jt^ r j^,rxSF{i)v,,{-i, 0) ® rySF{-i)vui, o)g,.(/). (iv.12) 



The color index contributions are already factored out in the above 

j^^)_^ = T^ir^,5^^®l = Cpl^l, (IV.13) 

j£^ = T^1Q1T^6^^ = -101 - ^101, (IV.14) 

41, = tH®IT^ = il0l - ^1®1, (IV.15) 

4 = T^lT^Ol = CplQl, (IV.16) 

4ti = TH®TH = il0l - ^1®1, (IV.17) 

jS = THqTH = h®l - ^101. (IV.18) 

Note that contributions should also be included where one or more of the gluon vertexes 
is replaced with V, from the clover term. We shall adopt an on-shell massless scheme in 
this section; quark mass and all external momenta are set to zero. 



A. Contribution from diagram (a) and (a') 

The one loop corrections Jj^y and Jj^y are the same as those to the bilinear (axial) vector 
current operator. Omitting the color factor the one loop contribution is given by 

l'^\ = Tv{V®A), 4^2 = Ta{V®A), (IV.19) 

I^^I = Ta{A^V), I^^;)=Tv{A^V), (IV.20) 

where V ® A = j^ ^ 7^75 and Ty/A is the one loop correction to the local (axial) vector 
current for the Wilson fermion 

TrT = j^ 7^^im(0, 1)Sf{1)VSf{1)VU-1. 0)^^.(0. (IV.21) 

Introducing the gluon mass A to the propagator G^y{l) in the loop we obtain [^ 

in the Feynman gauge, where /i2(r) is an integer given by 

h^{T) = 4(A), 4(\/), 16(P), 16(5), 0(T) (IV.23) 

for various Dirac channels. The finite constants Vr depend quadratically on the clover 
coefficients csw, and we write 

Vr = V^""^ + cswVr^'^ + c^w^r'^- (IV.24) 



The superscript (i = 0, 1, 2) means a correction of i-th order in csw, where i gauge interac- 
tions are replaced with that from the clover term. The numerical value of the finite part Vr 
has already been evaluated in Ref. [4] for various gauge actions and is given in Table |T] for 
Iwasaki gauge action. 

B. Contribution from diagram (b) and (b') 

The contributions from the diagrams (6) and (6') can be evaluated by using the Fierz 
rearrangement |5| in the spinor indices. 

Each one loop correction turned out to be the same as that to the (pseudo) scalar density 
and (axial) vector current operators given by (IIV.2ip . After carrying out the loop integral 
we obtain 

IvA+Av + ly'l+Av = -{Tv + Ta){VqA + AqV), (IV.25) 

IvA-Av + IvLav = -2 (Ts + Tp)iSQP-PQS), (IV.26) 

where the direct product © means 

(r r')„,,,, = (r)„, (r'),, (iv.27) 

and S Q P = 1 75. Performing the Fierz rearrangement again the vertex function is 
transformed into the same spinor structure as at the tree level 

IvA+Av + IvA+Av = {Tv + TA){y®A + A®V), (IV.28) 

IvA-Av + IvA-Av = {Ts + Tp){y®A-A®V). (IV.29) 

C. Contribution from diagram (c) and (c') 
Evaluation of the contributions from the diagrams (c), (c') is performed by using charge 







conjugation and Fierz rearrangement J5|. We use the representation of the charge conjuga- 
tion matrix C = 7270 and the relations 

CSMC~' = Spi-pf, (IV.30) 

CV,,{k,p)C-' = -V,,{p,k)^, (IV.31) 

CVl;\k,p)C-' = ~vi;\p, kf. (IV.32) 

9 



Using these relations and the Fierz rearrangement each quantum correction becomes the 
same as that to the (pseudo) scalar density and (axial) vector current operators 

I^v\+Av = 2Ts {SC-^ ® CP) - 2Tp {PC-^ ® CS) (IV.33) 

IvA-Av = Tv {VC-^ ® CA) + Ta {AC~^ ® CV) , (IV.34) 

^VA+Av = 2Tp {SC-' ® CP) - 2Ts {PC-' ® CS) (IV.35) 

IvA-Av = Ta {VC"' ® CA) + Ty {AC'' ® CV) . (IV.36) 

where the direct product ® means 

(r ® r')„,^,, = (r)^^ (r'),, (iv.37) 

Performing the Fierz rearrangement and the charge conjugation the vertex function is 
transformed into the same spinor structure as at the tree level without any mixing 

IvA+Av + ly'l+AV = -{Ts + Tp){V®A + A®V), (IV.38) 

IvA-Av + IvLav = -{Ty + TA){V®A-A®V). (1V.39) 

V. CONTRIBUTION FROM PENGUIN DIAGRAMS 

A. One loop correction 

In order to evaluate contributions from the penguin diagram we consider a four-quark 
operator of the following form 

Qfy = (T('=)),,^,,(rx®ry)„^^^, Yl ^^^Hsaad,^) {q,,q,s) , (V.l) 

q=u,d,s 



where k = 2n — l oi 2n and coefficients aq and T*^^) are defined in fllll.191) and fllV.3p . Then 



we evaluate the penguin diagram contribution to the amputated four quark vertex function 

h;XY = (^QxYSaa{Pl)dbl3{P2)qcy{P3)qdsiP4)/ > (V.2) 

which is given by the Feynman diagrams in Fig. |3l All the external momentum are set to 
in-coming direction and the internal gluon momentum is given hj p = pi + p2 = —{ps +^4)- 
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The vertex correction is given in the form 

(V.3) 

l2n;XY = Jpen Yl ^^^S^P, ^.) (r^ ® Vl,(p3, P4)) G^.Ip), (V.4) 

q=u,dyS 

where the color factor is given by 

'/pen=(-^l®l + ^l0l) (V.5) 

and we define the following loop integrals 



Iy^l\p,m,) = - I ^tr(ry5,(/-p)\/i,H+p, 05,(0), (V.6) 



" dH 

^^ (TxS.il - p)V^,{-l + p, l)S,il)Ty) . (V.7) 

We notice h-xv = 0. The contributions of lyL and Ixy u. should also be included where 
the gluon vertex in the loop is replaced with V-{^ from the clover term. 
We calculate the following loop integrals 

^yJ,M (P' ^i) = -^^ [lyl^iP. mg)) , /^^^(p, rrig) = -tr [i^id^ip, mq)) , (V.8) 
^yi^l(P> ^q) = IvI^iv^ "^9)7^^75, ^Iy;J,(p, "^<?) = Ivld^iv, mg)-^^, (V.9) 

I^{p,mg) = I ^^S^{l-p)V,,{-l+p,l)S^{l) (V.IO) 

according to the standard procedure of lattice perturbation theory |6|, i.e., we expand the 
functions in terms of the gluon momentum p^ and the quark mass m. 
The vertex functions satisfy the vector Ward-Takahashi identity |6| 

J]2sin^/;(p,m) = 0, (V.ll) 

5^2sin|i/;(^)(p,m)=0 (V.12) 

arising from the identity 

J22sm^V^,i-l+pJ) = Sp\l-p)-Sp\l), (V.13) 

J22sm^vl;\-l+pJ) = 0. (V.14) 

11 



Hence the loop corrections should be proportional to 

lim /J(p, m) oc {p'^6^^ - Pf,p^) . (V.15) 

a— >0 

We notice that a term proportional to a^^Pu is also allowed by the identity. However this 
term gives the same form of contribution as flV.lSP when expanded in terms of the external 



momentum and substituted into (IV. Sp and (IV. 9p . A detailed discussion shall be given later 
in Sec. IVITBI 

The only non-vanishing candidate is 



i2 



4iiyjp^-^) = ^w g(^Jg(^^^/ ";gy;,(p,^) 



+ C(a), (V.16) 

p=m=0 



which has logarithmic divergence and should be regularized with some infra-red regulator. 
We shall adopt the gluon momentum p^ as a regulator and the regularization term is defined 
by the similar loop integral 

^^,9{n^~l^)L',fl,Jl,pU, (V.17) 



where Ly/*^y (/,p)iR is the same integrand given in (IV.Sp . ( IV. 9^ but with all the Feynman 



'Y/XY;im\ 



rules replaced with that in the continuum and the quark mass set to zero. 
The loop integral is evaluated with a subtraction 

^y/Sy;M(P' "^) = ^Y/iY;M "') ~ 4/iY;M^n + ^Y/iYjP)^^- 



1 & 



2 



P(') (^ ^\ tP(') 



2^"^^ diap^)diap,) Vy/^y>^^^ - 'y/xY-AvU 



+ 0{a) 

p=m=0 



+4itYJp)m- (V.18) 



The first term is finite and can be evaluated numerically. 



a(ap.Map,) vC("'-^S(0).k) = ^ ((5.09290(43)) ^,„*„, 

- (1.88003(27)) (S^J.f + S^J^f) ) , (V.19) 



''' 'i:!>) ' iZW^k] = 0, (V.20) 



d{apa)d{ap,3) 



1 d^ 



Ci(0)-CJ(0)iK)7.75 



2 d{apa)d{api3) 

(V.21) 
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The second term has a logarithmic divergence and is calculated analytically 

(V.22) 



flU 



IvZ{p)iK = iI^UpW = ^C^(p)iR7.75. (V.24) 



rP(l)/ N 4? r. 2r \ 1 A ^^ 5\ p^ ^ 1 ^ 

h.„ (P)IR = Y^ [(-P V + PmP.) 3 ^1- ^ + 6 J - y ^- + ^^^ 

^2'i(p)m = 0, (V.23) 

We shall drop the last 1/a^ divergent term since the corresponding term is absent in the 
finite part since it is evaluated in terms of a derivative with external momentum. 

Contribution from the clover term is given by replacing the gluon interaction vertex with 
Vi^ in flV.6p and flV.7p . The loop integral has no IR divergence and can be evaluated 
numerically. 

-^S^(O) = TTT^f- (2.90088(27)) 5^,5,^ + (1.45031(13)) {5,J,p + 5,J,p)), 



d{ap^)d{app) "^-^ ' ' lOvr^ 

(V.25) 

iZui^) = 0, (V.26) 



d{ap^)d{appy^^''' 



d jP(2c)if.-. _ d jP(2c)if.^ _ 1 d tP(Ic), 



a(ap.)a(ap,)^^^^'-^°^ - 9(ap.)9(ap,)^^^i^°) " 2 9(ap.)9(ap,)^^--' ^^^^^^^^ ^^'^^^ 

B. Tree level contribution 

We consider the tree level contribution to the four quark operators given in Fig. HJ These 
diagrams may give a power subtraction with lower dimensional operators. We shall evaluate 
the amputated quark bilinear vertex function given by 



r(sub) _ /^(fc) 






irxY = {Q7ysac.{-P)d,,{p)) ^^^. (V.28) 

For each operator we have the following vertex function 

q=u,d,s 

+a'-;^6at4y\m,U + «i"^'5a6/??(m,)„^, (V.29) 

q=u,d,s 
+KafhabltY\^dU + A^c«i")<^afe4x'Vs)a/3, (V.30) 
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where we have two kinds of loop integrals in the above 

4t^M = j ^TxS,{l,m)Ty, (V.31) 

I?\m) = J^,tT{TySAl,m)) = (V.32) 

The latter vanishes for both operator Y = V,A. An explicit calculation gives 

4t)(m) = /(-'^)(am)75, (V.33) 

liv\m) = -I^^-''\am)j„ (V.34) 

J (271^1 +W{1, am)'' 
Substituting this result we obtain 

4n^A = -linHv = Na^:^Sa, (75)., (/^-^^ (m,) - J^^^) (m.)) , (V.37) 

which may be evaluated with an expansion in the quark mass 

a'' d[am) a 2a[am)'' Q d[amy 

(V.38) 



The numerical evaluation gives 



^ ./(sub)(o) = ^ (-21.46586(54)) , (V.39) 



d{am) 167r2 

5(^/"""(0) = ^(-14.9167(11)). (V.40) 

It may not be a good idea to expand in terms of the quark mass since these coefficients are 

rather large and furthermore d^ / d{amY I'^^"^\Q) term has an infra red divergence at m = 0. 

This contribution introduces a mixing with the lower dimensional bilinear operator (575^) 

multiplied with a mass difference {md — irts) as is given in (1111.32^ . It is clear from flV.35P 



that the mixing is due to the chiral symmetry breaking effect in the Wilson fermion. 

VI. RENORMALIZATION FACTOR IN MS SCHEME 

We renormalize the lattice bare operators Q\J^ to obtain the renormalized operator Q^L- 
We adopt the MS scheme with DRED or NDR. The renormalization of the operator is given 

14 



by 

Qg = zf,Q^ + zrQ^:: + zf^o^ ( vi. i) 

where Ql^l is the four quark operators on the lattice, Qjfcn is the QCD penguin operator and 
Ol^^ is a lower dimensional operator to be subtracted. Zf- comes from the gluon exchanging 
diagrams. Zf^^ is the contribution from the penguin diagrams. 

A. Gluon exchanging diagrams 

For gluon exchanging diagram we sum up all the contributions from three diagrams (a), 
(6), (c) and multiply by the color factor. Here we show its explicit form for the AS* = 1 
operators. 

I^vliiv^ = (^^^ {Ty + Ta) + ^ {Ts + Tp)^ (101) {V^A + A^V) 

+^ {Tv + Ta-Ts- Tp) (101) {V®A + A®V) (VI.2) 

4a1aT^ = (^^ (Ty + Ta) + ^ {Ts + Tp)^ (101) (V®A + A®V) 

+i (Ty + Ta-Ts- Tp) (101) {V®A + A®V) (VI.3) 

I^vl-'^v = [i^ {Ty + Ta)-^ {Ts + Tp)j (l0l) (V^A-A^V) 

+ ^ {-Ty -Ta + Ts + Tp) (101) (V^A-A^V) (VI.4) 

ll^X-'% = (^^ {Ts + Tp)^ (101) (V^A-A^V) (VI.5) 

From these vertex functions one can easily see that the one loop correction to the four quark 
operators is given in a form 

^one-loop ~ -^ij Qtrce^ (VI. Oj 

where QllL is a tree level operator. The correction factors are given by 

/\r2 _ 2 1 

11 - ^22 -^33 -^44 -^99 -^10,10" 2N ^^^ + ^^'' + 2iV 

g2 / ^2 _^ 2 o A^2 _ 2 1 \ 

^ ' ln(Aa)V— ^(W + K4) + — (V5 + V» , (VI.7) 



167r2 V N ' ' ' 2N ' ' ' 2N 

Tit = nf = ^{Ty + TA)-^{Ts + Tp) 
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„2 / A7-2 _ 4 N 1 

^, [-^^ In (Xaf + y (VV + V^) - — i+Vs + V^) ] , (VI.8) 



nf = nf = -1^^ (Ts + Tp) 

' -4— ^— ln(Aa)' + — -— (\/s + Vp)), (VI.9) 



167r2 V A^ ' ' 2N 



Tlat rrilat rrilat rrilat rrilat rrilat "^ /'rji , rji rji rji \ 

12 - ^21 -^34 -^43 - ^9,10 - ^10,9 - 2 '^^^ + ^^ ^ ^■^ ~ ^^'' 



- (6 In (Aa)' + Vv + Va-Vs- Vp) , (VI. 10) 



167r2 2 



^^ 1 



(-6 In {XaY -Vv~Va + Vs + Vp) , (VI.ll) 



16772 2 

where A is a gluon mass introduced for infrared regularization. 

The renormahzation factor is given by taking a ratio of quantum correction with that in 
the MS scheme multiphed with the quark wave function renormahzation factor Z2 

(zf^yu+T^) 

Zf.{^ia) = ^ ^-^ ^, (VI.12) 

ZfM) = T!t' - T., {1+3). (VI.13) 

The correction factor in the DRED MS scheme is given by 



mMS _ rpMS, _ rpMS> _ mMS _ mMS _ mMS _ / ^ " "^ ^ \ t/MS (\ti i /|\ 

-'ii — -'22 — -'33 — -'44 — -'99 ~ -'10,10 ~ 1 J^T ) ^ ! i^vi.i^; 



'Ar2 + 2^ 

T-iMS _ T-iMS _ rpUS> _ rpUS> _ TiMS _ rpUS> _ _OT/MS (\T-I 1 f-\ 

-'12 —-'21 —-'34 —-'43 ~ -'9,10 ~ -* 10,9 ~ "^^ ' l^Vi.iOJ 

Tjf = Tjf = (^^) V^^, (VI. 16) 

T^ = T^ = W^, (VI. 17) 

TW = TW = ^^^y^^ (VI.18) 

y-^J^(iog(g)+i). _ (VI.19) 

The quark wave function renormahzation factor is given by Ref. [4| and the result in the 
Feynman gauge is 

7MS\ Jl , 

'2 1 /,_\ _ 1 I y n ( 1^„/',,^^2 I vMS 



Z. 



lat 



{^a) = 1 + ^Cp (- log(/xa)2 + Sf _ Si) , (VI.20) 
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where 

Ef^(DRED) = --, (VI.21) 

El = Sf) + cswSS'^ + 4^^?. (VL22) 

The numerical value of S^ is given in table HTl 
Substituting the above results we have 

Zf^ifia) = Z^^ifia) = Zlsifia) = Zl^ifia) = Z^gifia) = Zfg io(/ia) 

Z|,i^^a) = Z',,ifia) = 1 + -|!^ (^-1 In i^^af + zl,^ , (VI.24) 

ZleiH = ZU^^a) = 1 + ^1^ ( ^^^^~^^ In {H" + 4^ , (VI.25) 

Zf2(/ia) = Zfi (/ia) = Zl^{^ia) = Zl^ifia) = ^|,io(/^a) = Zf^^gifia) 

^ {-^Inifxaf + zf^), (VI.26) 



Z|6(/ia) = Z^,,i^,a) = yI^ (3 In (fiaf + zl,) , (VI.27) 

Zl,ifia) = ZU^,a) = ^zl, = 0. (VI.28) 

4i = -^ ^^ {Vv + ^a) - ^ {Vs + Vp) + 2Cp (sr - Si) ,(VI.29) 

4, = ^^ -^(^y + ^^) + ^ (^^ + ^^) + 2Cf i^f" - Si) ' (vi-30) 

^6 = 4^^^ ^^ (V5 + Vp) + 2Cp (Sf - Si) , (VI.31) 

^fs = -3-^{Vv + Va-Vs- Vp) , (VL32) 

^6 = -42- (VI.33) 
The numerical result is given in table IIIII as an expansion in csw 

~9 _ ^ff(O) _|_ „ -9(1) I „2 ^9(2) f^j-T o^N 

for A^ = 3. 

We need to subtract the evanescent operators E^^^ in the MS scheme, which comes from 



QE 



the difference of dimensionality from four for gamma matrices in the operator vertex 
The evanescent operators in the DRED scheme is given by 
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Ei,3,9 = (1®1) (^li ® 7^ - lu ® lu) \, (VI.36) 

i?2,4,io = (101) (^7^ ® 7^ - lu ® 7.'') I (VI.37) 

E,j = (1®1) (^7^ ® 7^ - ll ® 1^^ I (VL38) 

^6,8 = (101) {^lu 7^ - li 7.'') ^, (VI.39) 

^(DRED) ^ ^(DRED) ^ ^(DRED) ^ _^;v^ (VI.40) 

p(DRED) _ p{DRED) _ p(DRED) _ 9 , . ,, .. , -. x 

-&22 - -^44 - -^10,10 - Y^^^' (Vi.41j 

p(DRED) _ p{DRED) _ p(DRED) _ q'^ /, .. .^x 

-^^12 --'^34 - -^^9,10 ~ 16^' (,Vi.4/:j 

p(DRED) _ p(DRED) _ p(DRED) _ q"^ ,, .. -oN 

^21 - ^43 - ^10,9 - ~Y^' (Vi.43j 

^(DRED) ^ ^(DRED) ^ _^2C^, (VI.44) 

E(r-)=E(r°) = ^i^, (VI.45) 

E(f -°) = E(°^'^°) = --|!^, (VI.46) 

where n is the dimension of the loop momentum. ^^^^ and 7,/ are n = (4 — e) and four 
dimensional gamma matrix with the chiral projection (1 =p 75). 
The conversion formula to the NDR scheme is as follows J9| 

(^?i)^°"=(^?i)°""°-^, (VI.47) 

(45)"°" = (4)°""° - ^, (VI.48) 

(-Ia)"°" = (-Ia)""""-^, (VI.49) 

(.?,)-- = (.?2)°""° - ^, (VI.50) 

(4a)"°" = (4a)°""° - I (VI.51) 

(4,)^D« = -3 (VI.52) 

with corresponding evanescent operators. The numerical value of zf, is given in t able HVl for 
NDR. 
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B. Penguin diagrams 

Taking a summation of the finite part and the IR divergent term the one loop correction 
from the penguin diagram is given by 

C-i;yA = c^^^Xtip) (-^1®1 + 101) (7m75 ® 1,) , (VI.53) 

C-Mv = «f r;t(p) (-^1®1 + 101) (7.75 ® 7.) , (VI.54) 

Cr^A = 0, (VI.55) 

Cr^v = I E 4"M ^S(P) f-4l®l + loO (7.75 ® 7m) , (VI.56) 



E < ) ^S(P) (-^1®1 + 101) (7.75 ® 7m) , 

^2 o 



Tltip) = if^^ (1^ «V + C) ' (VI-57) 

where we adopt the on-shell condition for external quarks 

- i^3 + mg = 0, i^i + mq = 0, p= -ps - p^. (VI.58) 

The finite part is expanded as 

^pe'n = Vi^l + cswV^ll (VL59) 

with coefficients given in table |Vl 

We notice that the above vertex corresponds to a four fermi operator of the form 

(Sa7M754) E ^'^bll^Qb), {Sal,,l5db) E ^'^bll^Qa) (VI.60) 

q=u,d,s q=u,d,s 

and is given by a linear combination of Q^^\ Q^'^\ Q^^\ Q^^\ which defines the penguin 
operator 

Q = [QvA+AV + QvA-AVj ~ jY [QvA+AV + QvA-AV ) ■ (VI. 61) 

The one loop correction from the penguin diagram to the four quark operators is written as 

Ql^ioo, = (Trr Qlt, (VI.62) 

where Qtree is ^^^ penguin operator at tree level. The correction factor is given by 

(Trr = il^^^ (1^ «v + o (VI.63) 
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with operator dependent factor 



C (g(^)) = 0, (VI.64) 

C (Q(2)) = 1, (VI.65) 

C (g(3)) = 2, (VI.66) 

C (Q(^)) = C (Q(«)) = Y. "f = ^/' (VI-67) 

q=u,d,s 

C (Q(^)) = C (Q(^)) = 0, (VI.68) 

c (g(«)) = c (g(i°)) = Y. ^f = ^« - T' ^^^-^^^ 

C {Q^^^) = -1. (VI.70) 
The correction factor in the MS scheme is given in a similar form 

<?ll-i„„„ = (Tr'^QSi, (VI.71) 

,^..p^_|LC(|!)(,„g)_5_^(^,,)J_ (V.,72) 

where the scheme dependent finite term is given by 



(NDR) ^g(2)^ ^ ^(NDR) (Q(2n-1)) ^ .^^ ^(NDR) ^g(2n)^) ^ q^ (yj 73) 

(DRED) ^g(2)^) ^ ^(DRED) (Q(2n-1)) ^ ^(DRED) ^g(2n)^) ^ 1 ^yj_74) 



for ra > 2. 

Combining these two contributions the renormahzation factor for the penguin operator 
is given by 



I r^TT ^ -i ^ ' 



167r2 3 



^r = -C-3-c(Q«)- (VI.76) 

Numerical value of the finite part is given in table I VII 

C. Subtraction of lower dimensional operator 

As was discussed in Sec. IVBI the lower dimensional operator 

OiTt' = -si^d (VI.77) 
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mixes with the four quark operators. The subtraction factor is given by 

ZS = -2af (/(^"^)(m,) - /(-^)K)) , (n = 3, 4), (VI.78) 

Z£"') = -2Naf (/(^'^'^Hm,) - &''''\m,)) , {n = 3,4), (VI.79) 

Zt-l = zt'^ = 0, (n = 1,2,5). (VI.80) 

We may be better to evaluate these factors nonperturbatively for numerical simulation. 

D. Mean field improvement 

The mean field improvement is given by subtracting the tadpole contribution in the 
renormalization factor and replacing it by a nonperturbative value u given in terms of the 
average plaquette u = P^'^ for example. The tadpole contribution resides only in Si of the 
quark wave function renormalization factor Z2. The mean field improvement works for the 
diagonal renormalization factor Zf^ from the gluon exchanging diagram. In the improved 
renormalization we shall use the renormalization factor 

u'zi'^^'^ (VI.81) 

instead of Zf^. Z^^ is given by replacing the finite term zf^ by z^^ in which the tad- 
pole contribution is subtracted. The csw dependent part is not affected by the mean field 
improvement. The numerical value is given in table IVlTI 

VII. 0{a) IMPROVEMENT COEFFICIENTS 

In order for the on-shell 0(a) improvement program to work one need to adopt the rotated 
field for the operator 

ar / 



V^c = [1 - y [z-i^D^ - (1 - z)mjj 7/>, (VII.l) 

^c = V' [1 - Y [-zi,D, - (1 - z)m)\ (VII.2) 
in addition to the improvement of the action. We shall set the on-shell condition 

{l^JD^, + m^ ^g = 0, ^^ (-7^5;, + m^) = (VII.3) 

for the quark fields and adopt 2; = for simplicity. The bare mass nig is defined by sub- 
tracting the additive mass correction from the bare Wilson fermion mass. 
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A typical form of the tree level improved four fermi operator is given by 

0iat = [I + — {1711 + 1712 + 1713 + m^)j (^iF^s) (i^s^'i^i) + 0{a^), (VII.4) 

which we shall adopt for our lattice bare operator. Each quark fields ipi has an incoming 
external momentum pi. In the following we set r = 1. 
The renormalization relation is given by 

Qg = z!,Q\ii + zrQ^- + zr'o!::^ 

with 0{a) subtractions, where Qiat is a tree level improved lattice bare operator for the 
K —^ TTTT decay. On.\a.t represents four fermi operators with wrong chirality given in (IVII.SP 
- fIVII.lOl) . i?'s are proportional to the quark mass rriq. Q's are dimension seven operators 
proportional to the quark external momentum p^. 

A. Contribution from gluon exchanging diagrams 

We shall evaluate the 0{a) correction for massive quarks in this subsection. The cor- 
rection has already been calculated in Ref. 3J] for gluon exchanging diagrams with massless 
quarks. 

For the gluon exchanging correction all the ten operators Q^^^ are not distinguishable but 
we have only four distinction 0|/J^_|_^y, where k takes even or odd for the color factor. So 
we shall evaluate the one loop correction to the following four operators 

Oi = o!;i^^^, 02 = o^\^^y, (VII.6) 

Os=0^\_^y, 0, = OP^_^y, (VII.7) 

for which we shall need six more operators to mix with at 0{g'^a) 

0, = 0fp_ps, 0, = 0fp_ps, (VII.8) 

07 = 0g,^p^, 0, = Ofp^ps^ (VII.9) 

O9 = Og, Oio = Og. (VII. 10) 

The fiavor structure shall take the form given in (lIII.2ip - (1111.25^ for a practical use in 
i^ — 7- TTTT decay. 
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We consider the following four fermi operator 

On = i^n) aabl3;c'yd5 (V'l;a,a^2;fe,/3) {^S^c^^i-As) , (VII.ll) 

(r---)„„...,. = (^'")*- i^y^^A-A^ VU„, . (VII.12) 

(^--4*.-,.. = (^'")-:- (S « P ± /> » SU^, . (VII.13) 

(4"=') =(r'")„s.,,(r»r) , (vii.m) 

V ^^ J aabl3;c^dS ^ ^""'^"^ V / a/3;7<5 

where T'^'^^ represents the color factor. 
The one loop contribution is written as 

^h,VA±AV = -^fc"'^''^ ( 1 + -a (mi + ma + ma + m4) j Iva±av (VII. 15) 

including the 0{g^a) terms, where J^' is a color factor given in Sec. IIV[ 

As was mentioned there the one loop correction is given in terms of that to the bilinear 
operator (IIV.21D 

(VII.16) 

4Z\v = {g^v'^ &A + VQ Gf ) + gT^ QV + Aq G^^'^) , (VII.17) 

4a1v = -2 {g^s"^ &P + SQ Gf^ - g'^^^ QS-PQ Gf )) , (VII.18) 
4aSv = -2((Gf )g-i ® GP) + {sG-'®GGf^) 

- (g'p'^^G-^ ® GS) - (PG-^®GGf^)\ (VII.19) 
4T']iv = {g^v'^C-' ® Ga) + (yG-i ® GGf^) 

+ (cf^G-^ ® Gf) + (aG-^ ® GGJ,^^^) , (VII.20) 

where Gp is a one loop correction to the bilinear vertex T with i-th and j-th flavor con- 
tributes for the internal quark line 

7^-^yif^{PiJ - Pi)SF{l - Pi,mi)TSF{l + Pj,mj)Vi„{-l - pj,pj)G^u{l). 

(VII.21) 

The quark mass and external momentum is kept non- vanishing here. 

The vertex correction can be expanded in terms of the quark mass and the external 
momentum up to 0{a) according to Ref. [lO| 

Gf!^'^ = TrV + T^™^ ^a (m^ + m,) T + T^''\a {p, + p,)^ F^, (VII.22) 
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where the vertex F^ for operator subtraction is given as r+ in Ref. 10|. We notice that 

0{g'^a\oga) term cancels by adopting the tree level improvement condition csw = 1 and the 

on-shell condition ( 1A.19P for the external quarks. All the coefficients are proportional to g^. 

Performing an explicit evaluation with csw = 1 the O(q^a) coefficients are given by 



<p) 



-Ca and T, 



(p) 



Cy, where Ca/v is defined in Ref. 



10( |. We notice that definition 



of the tensor operator cr^j,^ = ^ [7^^,71/] is different from that in Ref. 10|. Tp = Tg = 



.(p) 



for on-shell quarks. Tp is equivalent to Vp in Ref. (lO| and its explicit value can be 
reconstructed by 



Vr 



(1) 



-Br + bo 



T ---r^^o (VII.23) 

in the reference. We give numerical values in tables IVlnl IIXI for the Iwasaki gauge action in 
order to avoid confusion. The coefficients for the other gauge action can be reconstructed 
according to the example from the numerical tables in Ref. 10 1. 

Substituting into (1VII.16I) - (1VII.20I) the one loop correction is expanded as 



^(diagram) 



T(diagram)(g^) t (diagram) (g^ am) t (diagram) (g^ap) 

l\/A^Ai/ 'T -'■V A + AV + -'(, 



•■VA±AV — '^VA±AV '^ ■^VA±AV ~^ '^V A±AV ■ (VII. 24) 

Each term represents 0{g'^), 0{g^am) and 0{g^ap). 

In this subsection we adopt the lattice scheme (1IV.22I) with gluon mass regulator implicitly 
and evaluate the full vertex correction Gi, ■ ■ ■ , G4 for operators On=i^4. For example Gi is 
given by 



Gi 



1 + -a (mi + 1712+1713 + 1714 




(o) , Y- 

VA+AV ^ / J 

diagram=a,6,c 



T(diagram) 7-(diagram) 



(VII.25) 



with the tree level vertex ^ya+av 







According to Ref. p^ we rewrite the bare quark mass in terms of the renormalized one 



VTia = Zrr,m. 



mi'''qRi 



Z« = Gi.(-3L + z^), L 



167r2 



InA^a^. 



(VII.26) 
(VII.27) 
(VII.28) 



We then multiply the four quark vertex correction with the wave function renormalization 

1 + ^/4^) + \am,R (-1 + g^Z^^) , (VII.29) 



factor Z'J' Z'J' Z'J' Z'J' 

Vl V>2 i>3 i>4 



A/2 
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4'^ = Cf {-L + El) , (VII.30) 

4"^ = Cf (^L + ^Si ~z^ + J:^A . (VII.31) 



Here Si, 111 and Zm are constants introduced in Ref. 10| and is given in table|X]for csw 
The 0{a) and 0{g^ a log a) terms are canceled in 



,1/2^1/2^1/2^1/2^ _ ^si(lat) 



ZrzrZ^^Z^^G, = Zfp>T, + g'aB^.T, + g'aC.^, (VII.32) 



and one can extract the renormahzation factor Zf^^ and the 0{a) coefficient 5,^, B[^^ Cij. 
The renormahzation factor Zf ^ has aheady been used in Sec. I VI I to get that in the MS 
scheme. Although an explicit form of the 0{a) coefficient is given in the appendix [XJ we 
mention that the lattice bare operator Q\ll is given by 0„=i^4 and On,iat is given by 0„=5^io. 
The dimension seven operator Q^H is given by 

Q^'^ = (f,) ^ (?i;a,a^2;6,/3) (V^3;c,,V^4;d,^) (VII.33) 

V / aabpic-ydo 

in terms of the vertex given in flA.721) - (IA.79p . 



The same 0{a) coefficients appear in the renormahzation relation (JVII.SP for the MS 



scheme. These coefficients are written in terms of the one loop corrections Tp , Tp , 
El + S^ multiplied with quark masses and external momentum as will be given in the 
appendix. 



B. Contribution from penguin diagrams 

The one loop correction from the penguin diagram to the improved operator (1VII.4I) is 
given by a slight modification of the one loop vertex (JV.3P and (IV.4P multiplied with the 
tree level improvement factor 

M + - {nid + mj j (1 + amg) l2n-i/2n;XY- (VII.34) 

However this factor shall be canceled with the tree level contribution of the wave function 
renormahzation factor flVII.29p and we abbreviate it. We shall evaluate fIV.Sp . fIV.Qp by an 
expansion in the quark mass and the external momentum. 

Before performing the expansion we make use of the Ward-Takahashi identity (JV.lip . 



(IV.12p . We notice that the identity is valid on the lattice with a non- vanishing quark mass 
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and external momentum and impose a restriction on the one loop correction f IV.lOp 



^/fb. f^g) = (p^^mA - P^^Px) Ixip, rrig) + a^xpxl^ip, mg), (VII.35) 

p = 2sin-. (VIL36) 

Taking into account the Lorentz covariance each term is expanded to give the 0{a) 
contribution 



^xiP, mg) = l^lx + mql^-ix + Pxlps + Pa(ya\lpA^ (VII.37) 

"q^m I ""q'-mm'^ Potla^-p ~T~ P ''pp- 



Z^(p, m,) = /^ + rrigC + ^ICm + Pc.lX + pXv (VII.38) 



We substitute this expansion into a typical correction term (1VII.35|) and substitute further 



into the one loop penguin contributions flV.81) and flV.9l) . After a short algebraic calculation 
we find that only three terms contribute to the one loop penguin diagram 

Ivl^l (P' '^'i) = -4 ( Vp' - P,Pu) {l"" + ? + arUgC) , (VII.39) 

C;!(P,'^.)=0, (VIL40) 

4aUp^ ^g) = IavUp^ "^9) = 2 ^5,M (P' ^ih-^^^- (VII.41) 
From the clover term contribution 

j;W(p,m,) = J -^^SAl-p,mg)v};\~l + p,l)Sp{l,mg) (VII.42) 
we have the similar form of correction 

Iv!^^ iP, rn,) = -4 {S,.p' - p.p.) csw (/^(^^ + ?(^^ + anigC^^) , (VII.43) 

^vS (P' "^9) = ^™ (P' "^9) = 2 ^vJl'^ ^^' "^9)7^^75- (VII.45) 

The 0{g) terms have already been evaluated in Sec. |V] and we put the same result for 
the notation used here 

/f = (/^ +r;) = ^^ (-InaV + 1.7128269(84)) , (VII.46) 

im ^ (^iPio ^JPio^j ^ _^1 (1.087821(3)) . (VII.47) 

We notice that both the 0{ga7nq) coefficient l^ and Im has a logarithmic IR divergence. 
The same regularization scheme is also used here as was adopted in Sec. |V]and the coefficient 
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is given by 

C = Y^^ (^^«V - 3.59121(36)) , (VII.48) 

e^ = Y^^ (-^ InaV + 0.74846(23)) . (VII.49) 

We substitute these results into the penguin contribution flV.3l) and flV.4l) including all 
the contributions up to 0{g'^a) 

hn-iyA = l2n~l;AV 

= Jpcn Yl <^ {^vaUp^ ^«) + ^VA% rn,)) ® (v^^.P^) + vt\ps,p,) 

q=d,s 

xGf,u{p) 



= 2«(7 Jpen Y, "S"^ ( Or + csw/f ^'^) 7^.75 ® 7/. 

q=d,s 

+m, {li + csw (/f + C^'O + 4w^r^'^) 7m75 ® 7/. 

+ (1 - csw) (/f + csw/f^'^) 7m75 ® ^ (P3 - ^4)^) , (VII.50) 

hn-yA = 0, (VII.51) 

/2n;AV' = ^pen Y «i"^ (^^ ^P' "^<?) + ^^^^P' "^9)) 

X (^7^75 ® Vlp(p3,P4) + 7^75 ® Vl^\p3,p^)j G^p{p) 
= AigJ^,^ Y "g"^ ( {}i + csw/f^'^) 7m75 ® 7/. 

q=u,d,s 

+m, (l^ + Csw (^f + C"^) + ciw^f^'^) 7m75 ® 7/. 

+ (1 - Csw) (/f + csw/f^'^) (7^75 ® ^ (P3 - ^4)^) \ ■ (VII.52) 

Here we made use of on-shell conditions for the external momentum 

^(^l+^2)=0, ^(^3+^4)=0, (VII.53) 

{P3-P4)^P^ = 0, (VII.54) 

i(^f^u (ps + Pi)„ = i{P3- Pi)^ - 2m57^. (VII.55) 



We notice that the mixing with an operator 7^75 (8>i (ps — P4) drops if we set the improve- 
ment coefficient csw = 1- The 0{g'^a\oga) term in /^ and Im cancels in a combination 
^m + Csw ( ^f + Im ) for Csw = 1- The 0{g'^a) improvement is accomplished just by shifting 
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l^pen in dVL59|) as 

Vlt ^ Vlt + m, (0.04210(43)) . (VII.56) 

VIII. CONCLUSION 

In this paper we have calculated the one-loop contributions to the renormalization factors 
for the parity odd four quark operators, which contribute to the K ^ rnr decay amplitude, 
for the improved Wilson fermion action with the clover term and the Iwasaki gauge action. 
The operators are multiplicatively renormalizable without any mixing with wrong operators 
that have different chiral structures except for the lower dimensional operator. The 0{g^a) 
improvement coefficients are also calculated for massive quarks imposing csw = 1 and the 
on-shell condition. 
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the one loop correction is 

Appendix A: 0{a) contribution from gluon exchanging diagram 

In this appendix we evaluate the one loop correction from the gluon exchanging diagrams 
with non- vanishing quark mass and the external momentum. 

We start from the one loop correction flVII.16|l - flVII.20|l and substitute into flVII.22|) . The 



correction is expanded as (JVII.241) . Explicit form of the 0{g^am) terms are given by 

1 

"2 



Wa±av 



j-{b){g^am) 
^VA+AV 



±1 (^^rn, + m,) Tf ^ + {m, + m,) T^")) (A ® \/) , (A.l) 

_1 ((mi + m,) Tir^ + (m2 + m^) T^^) {V A) 

A [{rn^ + rni) Tf^ + (m^ + m^) T^^) {AqV) , (A.2) 
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'■VA+AV 



Tic){g am) 
^VA-AV 



+ ((mi + m^) Tp^ + (m2 + m^) T^'")) (P 5) , (A.3) 
/W(^'r^ = ({m, + ma) Tf ^ + (m^ + m4) T^T^) (5^-^ ® CP) 

- ((mi + mg) T^"^ + (m2 + m4) 4"^^) (PC-^ ® CS) , (A.4) 

^ ((mi + ma) T^ ^ + {m, + m,) TJ™)) (VC-^ ® CA) 

+i ((^^ + ^3) TJ™) + (^2 + m^) 4™)) (AC-i ® C\/) . (A.5) 

Using the Fierz transformation with or without the charge conjugation we rearrange the 
spinor structure 

{yQA) = -]-{y®A + A®V)-S®P + P®S, (A.6) 

{^AQV) = -]-{y®A + A®V) + S®P-P®S, (A.7) 

{SQP) = ^[s®P + P®S-V®A + A(E)V-f®T^, (A.8) 

{PQS) = ^[s®P + P®S + V®A-A®V-f®TY (A.9) 

{SC-^ ® CP) = j{s®P + P®S-V®A~A®V + f®TY (A.IO) 

{PC-^ ®CS) = ^[s®P + P®S^V®A + A®V + f®TY (A.ll) 

(\/C-^ ® CA) = -- (V ® A- A®V) - S ® P + P ® S, (A.12) 

(AC-i ® CV) = -- {V (E) A- A(E)V) + S ® P - P ® S. (A.13) 

Here we notice there appears an operator mixing with wrong chirahty proportional to the 
quark mass difference and chiral symmetry breaking effect {Ty — Ta) or (T5 — Tp). 
The 0{g^ap) terms are given by 

4A&y = 4'^^« (Ps + Pd, (7m ® 75) ± T^I^ia (pi + P2), (75 ® 7^) 

+Tl,''\a (pi + P2), (^M. ® 7m75) ± T^^'^^a {ps + Pa), (7^75 ® ^m^) , (A.14) 
4'S? = -Tf^a {p, + P3), (7. 75) - T^I^^a (pi + p^), (75 7.) 

-T^\a {pi + Pi), {a^, 7^75) - T^^\a {p2 + ps), (7/^75 (^t^u) , (A.f 5) 
IvA-Iv = -2Tf ^o iPi + Pa), {l, 75) - ^T^ia {p, + ^3), (1 7^75) 

+2T'^j^\a (pi + p4)^ (7^75 1) + '^Tf^a {p2 + Ps)^ (75 7^) , (A.16) 

29 



41?;? = 2T^^\a {p, + P3), {l,C-' ® C75) + 2Tl^\a fe + P4), {iC^' ® ^7,75) 

-24^)^0 (pi + ps)^ {l,l5C'' ® CI) - 2T^''\a (p2 + P4)^ (75^-^ ® C7m) , 

(A.17) 
4K? = Tj^^a {p, + p4), {l,C'' ® C75) + T^^^a {p, + ^3), (75^^'^ ® C7,) 

+4^^^a (pi + P3), {(r,uC-' ® C7^75) + T^'^^o (P2 + P4), (7^75^"^ ® Ca,,) . 

(A.18) 

Using the momentum conservation relation Pi + P2 + Ps + P4 = and the on-shell condition 

{i:^i + rrii) ipi{pi) = 0, ^i{pi){-i:^i + mi) =0, (A.19) 

{z]^i + m,) C-^V^f (p,) = 0, ^Pfip^)C {-i^, + m,) = (A.20) 

we rewrite the correction 

^A&y = -T^A^ (^1 - ^2) (5 ® P) T Tj^a (m3 - m,) {P ® 5) 

+T^^^ia (pi + P2), (a^. ® 7^75) ± Tl.^'^a (p3 + P4), (7^75 » f^^^^) , (A.21) 

IvA+lv = T^A<^ (^1 - ^4) (5 P) + Ti^^a (-m2 + m^) {P Q S) , 

-r'^ha (pi + pi)^ {a^y 7^75) - T^\a (p2 + P3)^ (7^75 (^tiv) , (A.22) 

/W(^';^) = -2 (rf a (mi - m,) + T^f^a (m2 + m,)) {S P) 

+2 (^Tjf^a (mi + m4) + T^i^ {-m^ + m3)) (P 5) (A.23) 

/(f)^(s^;^) ^ 2 (t^ a (mi - m3) + Tf a {m^ + m,)) {SC'^ ® CP) 

-2 fe^a (mi + m3) + T^^^a {-m^ + m4)) (PC'^ ® CS) , {AM) 

jW(^^;^) = -Tj^^a (mi - m3) {SC-' ® CP) - T^^a (-m2 + m4) {PC-' ® CS) 

+T^^\a (pi + p3), {a,,C-' ® ^7^75) + 4"^2a (p2 + P4), (7^75^^"' ® Ca,,) . 

(A.25) 

Taking summation of all the contributions for VA + AV and VA — AV and multiplying 
the wave function renormalization factor we extract the 0{g'^am) coefficients in the 

Pii = g\M ^C, (Si + S«) + ^C, - ^) (t^ + Tf) ) + ^ (t^^^ + T^^')^ 

+g^^ (T^^^am(i4) + tI^^gm'^ , (A.26) 

B,, = g\M]^ (rlr^ + T^r^ - Tf ) - Tir^) - g^ (t^ am(i4) + T^j^^uM) , (A.27) 
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B[, 



B'lQ 



5,3 = -9'am^v2) {c, (t^^^ - T^^) + ^T^^) , (A.28) 

5i4 = (7'am(i2)iTt\ (A.29) 

g'am,, (-^ (t^ ^ - Tf )) + 2C^Ti^)) , (A.30) 

/ami4i(Tf)-Tf)), (A.31) 

+ {^Cf + ^^ Tt^am(i3) + ^ (T^"^am(i2) - T^f ^am^a)) ) , (A.32) 

B'i9 = 9^ ( ""^(13)^ (^i"^ - ^p"j - l^Tfamn^i) + — (^T^^^am(i2) - T^^^am(i3) j j , 

(A.34) 
i?;,io = / (-«^(i3)^ (tI'"^ - Tp^) + ^Ti^^«^(i3) - (Ti^^am(i2) - T^^ am^,,))^ , (A.35) 

^21 = g^aM- (t^^ + TJ'") - T^ ^ - TJ,")) - g^ (T^^^am(i4) + T^^^aikf) , (A.36) 

B22 = g'aM {c, (S, + S«) + (c^ - ^) (Tf ) + Tf) ) + ^ (t^™) + TiT^)) 

+g^^ (T^^)am(i4) + T^^^cM) , (A.37) 

^23 = -^/m(i2) (t(-) - Tf )) , (A.38) 

i?24 = ^/^(i2) (Tf ^ - Tf)) + g'am^v2)CpTf, (A.39) 

5^5 = g'am^i,)T^/\ (A.40) 

5^6 = /CpmM (Tf ) - Tf )) - /am(i4)^Tf , (A.41) 

B',, = -/am(,3)i (t^-) - Tir)) + / (t^ am(i3) - T^ am(i2) + Tif)am(i3)) , (A.42) 

+/ (- (c^ + ^) Tf am(i3) + ^ (Tf am(i2) - Tif)am(i3))) , (A.43) 

B',, = -/am(i3)i (tI™) - Tir)) - / (T(^)am(,2) - TJ^^am^.s)) , (A.44) 

i?2,io = /a^(i3)^ (tI'") - T^ ) + / (c^Tt)am(i3) + ^ (Tl^)am(i2) - Ti^^am^.s))^ , 



(A.45) 
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-B31 

= —n' 

2 



/am(,2)C^ (Tf ) - Tf )) + a'^Tfam^,,^, (A.46) 



i?32 = -(7'^Tt^am(i2), (A.47) 



-/I (Tf am(i3) + T^f^aM) , (A.48) 

534 = -g'aM]^ (rf ) + Tf ) - Tf ) - Tir^) + g" (t^ am(i3) + TJf^oM) , (A.49) 

B',, = -g'am^u)^ [rlr^ - Tf )) - /2C^Tt)am(i3), (A.50) 

i?^6 = ^^a^(i3)^(Tf^-Tr), (A.51) 

-/ ((2^^ - ^^ Ti^^am(M) + ^ (T^"^am(i4) - T^'^amd^))) , (A.52) 

^38 = -fl'^«"^(i4)2 (^i™^ - T^'j + 5'^ ( -2^A ^""^(14) + T^^^am(i4) - T<l^^am(i2) j ,(A.53) 



7-)/ 2 -'- frri'm) rr(in) 



(A.54) 
(A.55) 



^3,10 = fi'^«"^(i4)2 (^5"^ - TP) - 9'^ ( l^Tf'am^u) + (^T^*'^am(i4) - T^^^am(^u) 

541 = -/am(,2)i (Tf ) - Tf )) - /^Tf am^^), (A.56) 

i?42 = g'am^u)^ (t^ - Tf) ) + g^^T^^^am^,,,, (A.57) 

543 = 0, (A.58) 

544 = g'aMCp ((Si + si')) + (rf ) + 4™))) + /2C^ (t^ am(i3) + TJ^^aM^ , (A.59) 

i?;^ = /am(i3)i (t(J") - Tf) ) - /Tf am(i3), (A.60) 

5^6 = -/«^(i3)^ (^r - ^f + a'^^A^am^.s), (A.61) 

5^8 = -/am(,4)C^(Ti™^-T(r)) 

+^' (^Tt^am(i4) + 2Ci. (Tjf )am(M) - T^""^ am^,^)) ) , (A.63) 

5^9 = -(72llTt)am(i4), (A.64) 



32 



T~>l 2 r^ I rriirn) rriU 

^4,10 = 9 am(u)CF [Ts - ^p 

+g^ (^T^fam^u) - 2Cf (rj^^am^u) - T^^^amda))") . (A.65) 

where the g^ dependence is shown exphcitly. The quark masses used here is given by 
1 



M 



(mi + r?7,2 + T^s + ^4^ 



R ' 



™iij) = -^i^i + ^j -^i' -^j')r^ {i'J'} = {1,2,3,4}- {i,j}. 



(A.66) 
(A.67) 



The 0{g^ap) coefficients are given as 



Cii — —g CpTy , Ci2 — 0, Ci3 — —g iz^Ty , Cu — g -Ty , 



.(p) 



G 



21 



-.2'^rT^ip) r^__ _ „2 1 rp{p) Q 



g^i-T^'>, C,, = g'—T, 



2N 



■V ' 



23 



0, C24 — g CpTy , 



.(p) 



C'SS — g CpTy , C36 — 0, C; 
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-g Jpj^v ' Css - g -Ty , 



(A.68) 
(A.69) 
(A.70) 



,1. 



<p) 



C45 - g ^Ty , C46 - -g 2^- 

with mixing operator vertex defined as 



:T, 



(p) 



y , C47 — 5? -Ty , 



2lTn(p) ^_ _ ^2 1 Tn(p) 



48 = -^'^^r (A.71) 



f 1 = 


= l^lia{pi +P2)^ 


y^^u ® If^lb) - (7a. 75 ® (y^^u)) , 






(A.72) 


f 2 = 


= lQUa{pi+p2), 


i(^,,u ® 7^75) - (7m75 ® (rf,^)j , 






(A.73) 


f 3 = 


= l^lia{pi +P4)^ 


{(yf^u 7^75) - (7^.75 f^Mi.) j , 






(A.74) 


f 4 = 


= lQlia{pi+pi)^ 


{(^f^u 7/.75) - (7m75 CTf^u)) , 






(A.75) 


rs = 


= l^lia{pi +P2)^ 


{(Tf,^ 7^75) + (7m75 cr^^) ) , 






(A.76) 


f 6 = 


= lQlia{pi +P2)^ 


(f^Mi. 7^75) + (7^75 (^f^u) ) , 






(A.77) 


Fr = 


= l^lia{pi +P3)^ 




'®Ca^.))' 


(A.78) 


f 8 = 


= lQlia{pi+ps)^ 




'^®Ca^O)- 


(A.79) 
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TABLE I: Finite part Vr for bilinear operators [4]. Coefficients of the term c^T^in = 0,1,2) are 
given in the column marked as (n). Terms proportional to Cg-^^ are zero for pseudoscalar P. 



(0) 


V 

(1) (2) 


A 

(0) (1) 


(2) 


S 
(0) (1) 


(2) 


P 

(0) (2) 


(0) 


T 

(1) 


(2) 


6.275 


-1.725 0.637 


3.367 1.725 


-0.637 


2.533 6.902 


-0.293 


8.348 2.254 


4.615 


-1.150 


-0.327 



TABLE II: Finite constants for quark self-energy [J]. Coefficients of the term Cg^(n = 0, 1,2) are 
given in the column marked as (n). Tadpole contribution is also listed. 



(0) tad (1) (2) 



4.825 7.482 -1.601 -0.973 



TABLE III: Finite part zf- of the renormalization factor from gluon exchanging diagrams. The 
DRED scheme is adopted. The color factor is set to A^ = 3. Coefficients of the term Cg^(n = 0, 1, 2) 
are given in the column marked as (n). 



4i 

(0) (1) (2) 


(0) (1) (2) 


(0) 


(1) 


(2) 


^9 

(0) (1) 


(2) 


-23.596 3.119 2.268 


-25.183 5.420 2.922 


-18.041 


-4.933 


-0.020 


-2.381 0.451 


-2.020 
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TABLE IV: Finite part zf- of the renornialization factor from gluon exchanging diagrams in the 
NDR scheme. The color factor is set to A^ = 3. cgw dependent terms are the same as that in the 
DRED scheme (n = 1,2). 



^9(0) 


^9(0) 


^9(0) 


^9(0) 


^9(0) 


^9(0) 


-24.096 


-25.350 


-19.708 


-4.881 


-1.120 


-3 



TABLE V: Finite part in the penguin diagram contribution on the lattice. Coefficients of the term 
Cg-^(n = 0, 1) are given in the column marked as (n). 




TABLE VI: Finite part in the renornialization factor for the penguin diagram contribution. Coef- 
ficients of the term c^Y^{k = 0, 1) are given in the column marked as {k). 



zr"(DRED)(0) 


zr'(NDR)(0) zP^:^i(NDR)(0) zP7(NDR)(o) 


(^Pen)(l) 


-0.2039 


1.0462 1.0462 0.0461 


1.0878 




FIG. 1: One gluon interaction vertex, k and p represent incoming momentum into the vertex. 
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TABLE VII: Finite part zf^ of the renormalization factor from gluon exchanging diagrams for 

the mean field improvement given with the tadpole subtraction in the DRED and NDR scheme. 



scheme 


g(MF)(0) 
^11 


g(MF){0) 
^55 


9{MF){0) 
^66 


DRED 
NDR 


-3.644 
-4.144 


-5.231 
-5.398 


1.911 

0.244 



TABLE VIII: Coefficients of the 0{g'^ap) correction to bilinear operators for the Iwasaki gauge 
action. The improvement coefficient is set to its tree level value csw = 1- 



;i67r2)Tj') (167r2)T^P) {W7t^)T^;'^ {16tt^)T^''^ 



0.4519(22) -1.1478(67) 



TABLE IX: Coefficients of the 0{g^am) correction to bilinear operators for the Iwasaki gauge 
action. The improvement coefficient is set to its tree level value csw = 1- 



167r2)ri") (167r2)T,i™^ {m7:^)T^r^ (levr^lT^) 



)^A 



^V 



-0.9764(47) -0.995(13) -1.17592(48) -4.1201(36) 



TABLE X: Coefficients of the 0{g'^) and 0{g^a) correction to the quark propagator for the Iwasaki 
gauge action. The improvement coefficient is set to its tree level value csw = 1- 



gauge action ci 


(167r2)Si (167r2)sf) {IQt^^)z^ 


-0.331 


2.25022(44) -9.9266(69) -11.395(53) 
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a,a 



/3,b 



7,c 





b' 



FIG. 2: One- loop vertex corrections for the four-quark operator (gluon exchanging diagram). 
a, (5, 7, 5 and a, 6, c, d label Dirac and color indices respectively. 
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I — p 



l-p 




I — p 



I — p 




FIG. 3: One-loop vertex corrections for the four-quark operator (penguin diagram), a, (3,^,5 
and a,b,c,d label Dirac and color indices respectively. All external momentum pj's are in-coming 
direction. 




t 



q2 




7T<5 



FIG. 4: Tree level contribution to the AS" = 1 four-quark operator, a, /3 and a, b label Dirac and 
color indices respectively. An external momentum p is in-coming direction 
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